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This  report  investigates  curve  fitting  for  the  empirical  explanation  of  shower  wastewater  treatment 
by  ultrafiltration.  A two-part  model  is  proposed  where,  initially,  the  system  is  at  unsteady  state 
caused  by  the  formation  of  the  boundary  layer.  l,ater,  the  system  comes  to  a steady  state  when 
the  rates  of  boundary  layer  formation  and  decay  reach  constant  values.  The  unsteady  state  portion 
is  described  by  an  exponential  equation  and  the  steady  state  portion  is  described  by  a linear  equa- 
tion. 
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Block  20  (continued): 


The  model  was  evaluated  on  two  separate  ultrafiltration  systems  with  each  system  utilizing  a dif- 
ferent membrane  configuration.  These  configurations  were  the  spiral-wound  configuration  which 
fit  the  model  as  stated  and  the  hollow-fiber  configuration  which  fit  only  the  unsteady  state  portion 
of  the  model  due  to  daily  removal  of  the  boundary  layer. 

The  study  concludes  that: 

a.  A single,  empirical  equation  cannot  describe  permeate  production  by  ultrafiltration  of 
shower  wastewater. 

b.  A model  based  on  unsteady  and  steady  state  phenomena  of  concentration  polarization 
can  empirically  describe  the  ultrafiltration  process. 

c.  The  hollow-fiber  system  can  be  described  by  a single,  empirical  equation  of  exponential 
order  because  the  system  as  studied  was  never  allowed  to  reach  steady  state.  Therefore,  the  un- 
steady state  equation  described  the  system’s  operation. 

d.  The  spiral-wound  membrane  system  can  be  described  empirically  by  an  exponential 
equation  during  the  boundary  layer  formation  and  by  a linear  equation  during  steady  state  opera- 
tion. ■ 
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PROCESS  DESIGN  FOR  TREATING 


SHOWER  WASTEWATER  BY  ULTRAFILTRATION 


I.  INTRODUCTION 

1.  Previous  Investigation.  Since  the  first  description  of  the  osmotic  pressure 
phenomenon  by  Abbe  Nollet  in  1748,  scientists  and  engineers  have  been  intrigued  by 
the  properties  of  semipermeable  membranes.*  In  the  past  few  decades,  membrane  de- 
velopment has  made  membrane  separation  processes  practical.^  As  might  be  expected, 
membrane  processes,  particularly  the  ultrafiltration  process  which  is  shown  in  Figure  1 . 
are  playing  an  increasing  role  in  water  and  wastewater  treatment.  The  ultrafiltration 
process  depends  upon  a pressure  driving  force  and  a membrane  permeable  to  a com- 
ponent (B)  in  a lic]uid  solution  and  impermeable  to  the  other  component  (A)  in  the 
solution  (Figure  1 ). 


SOLUTION  OR  SUSPENSION  OF  A AND  6 UNDER  PRESSURE 


• • • 

i:  • 


• • • 


• • 


^ RETENTATE 
^(CONCENTRATED  A) 


MEMBRANE 


• • 


PERMEATE  (B) 


Figure  I.  Ultrafiltration  process. 

When  used  for  wastewater  treatment,  membranes  are  chosen  that  will  retain 
impurities  to  be  removed  from  the  waste  stream  on  one  side  of  the  membrane  and  col- 
lect water  of  desired  quality  from  the  other  side.  Therefore,  the  quality  of  effiiient  is 


* I'.  Daniels  and  R.  A.  Albcrty.  I^hysiail  ('.hi^misiry.  3rd.  t-d.,  Jolin  Wiley  and  Sons,  Inc.,  New  York,  1966. 

^ H.  C.  Gooding.  “A  Theoretical  and  lAperimental  Comparison  of  Unsteady  and  Steady  State  LTttat'iltration  in 
Thin  Horizontal  Channels."  M.  S.  Thesis,  CTemson  University.  1972. 


fixed  as  a function  of  tlie  membrane  structure  rather  than  as  a function  of  tiie  system 
operational  parameters. 


The  most  significant  variable  in  tlie  ultrafiltration  process  is  water  produc- 
tion expressed  as  fiux  with  the  units  of  fiow  (volume)  per  unit  membrane  area  per  unit 
time. 


System  optimization,  therefore,  is  done  by  maximizing  flux  and  minimizing 
the  rate  of  decline  of  fiux.  This  fiux  decline  is  due,  in  part,  to  a concentration  buildup 
of  the  impermeable  portion  of  the  bulk  solution  at  the  membrane  surface  (concentra- 
tion polarization).  The  system  is  configured  as  shown  in  figure  2 with  the  bulk  solu- 
tion velocity  parallel  to  the  membrane  surface,  thus,  reducing  this  buildup. 
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Figure  2.  Configuration  of  an  ultrafiltration  membrane  system. 

The  purpose  of  this  investigation  is  to  develop  a mathematical  model  that  de- 
scribes the  treatment  of  shower  wastewater  by  ultrafiltration.  The  model  is  developed 
using  transport  principles  for  a general  ultrafiltration  system  and  is  then  used  to  evalu- 
ate actual  ultrafiltration-system  operation. 

II.  LITI  RATURI-;  RHVII  W 

The  two  membrane-separation  processes  that  use  a pressure  gradient  across  a 
membrane  arc  reverse  osmosis  and  ultrafiltration.  They  difler  in  that  the  solute  and 
solvent  molecular  dimensions  in  reverse  osmosis  are  generally  of  the  same  order  ol  mag- 
nitude; whereas,  in  ultrafiltration  the  solute  molecular  diameters  are  at  least  ten  times 
greater  than  the  solvent  molecular  size. 
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2.  C onci’iitration  i'olari/aliun.  In  a continuous  operation  ol  reverse  osmosis  or 
ultratiltration  systems,  a phenomenon  commonly  relerred  to  as  concentration  ]iolari/a- 
tion  often  represents  tlie  rate  major/rate  controllin;:  step.^  ( oncenlration  polarization 
is  tlescribed  as  tlf  accumulation  of  nontransferred  solute  in  the  houiular>'  la>erol  the 
membrane  as  shown  in  I'igure  3 This  phenomenon  is  the  result  ol  the  concentration 
of  the  solute  at  the  membrane  surface.  As  the  solute  concentrates,  a concentration 
gradient  is  established  resulting  in  a back  diffusion  of  the  solutes  aiul  reducing  the 
transport  of  the  solvent. 


MEMBRANE 

MEMBRANE  BOUNDARY  LAYER 


Figure  3.  Causes  of  concentration  polarization. 


Solute  may  be  transported  away  from  the  membrane  surface  by  molecular 
diffusion  perpendicular  to  the  membrane  (lateral  direction)  and  by  forced  convection 
of  the  stream  How  parallel  to  the  membrane  (axial  direction),  liventually,  steady  state 

^ N.  l,dk.shiminarayanaiah.  "Transport  Phenomena  in  Artificial  Membranes.” f.'liem.  /terieus.  6.S.  No.  5.  49l-5h.*i 
(1965). 
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IS  achieved  , and  the  rate  of  Uansport  of  solute  to  tl\e  membrane  surface  becomes  equal 
to  ifie  rate  of  transport  away  from  the  surface  if  the  process  input  variables  remain 
imdiangccl. 

( onct'iitratioii  polarization  can  have  two  detrimental  ellects  on  the  mass- 
transfer  rate  in  an  ultrafiltration  or  a reverse  osmosis  system.  The  lirst  detrimental  et- 
fect  is  seen  by  considering  the  familiar  eciuation  relating  mass-transfer  rate,  driving 
force,  and  re.sistance: 

Rate  = (Driving  l-'orcef/f  Resistance),  (1) 

or.  more  specifically : 

MTR  = (DP-1>1M)/R  <-• 

where: 

MTR  = mass-transfer  rate 

DP  = transmembrane  pressure  drop 

Pl\l  = osmotic  pressure  of  nontransferred  solute  at  the  membrane 

surface 

R = resistance  to  mass  transfer. 

Concentration  polarization  reduces  the  effective  driving  force  (DP  - PIM)  by  increasing 
the  concentration  of  tiie  nontransferred  solute  in  the  boundary  layer  and  thus  increas- 
ing the  osmotic  pressure. 

The  second  detrimental  effect  is  caused  by  the  high  solute  concentration  at 
the  membrane  surface  precipitating  and  lorming  a gel  layer.  The  thickness  of  this  gel 
layer  will  increase  until  mass  transfer  through  the  gel  becomes  the  rate-limiting  step. 

In  any  event,  once  concentration  polarization  is  controlling,  increasing  the 
pressure,  or  driving  force,  on  the  system  will  not  increase  the  mass-transfer  rate.  1 he 
higher  pressure  causes  a thicker  layer  of  gel  and,  hence,  a greater  resistance  to  flow. 
The  goal  in  process  design  is.  therefore,  to  minimize  the  thickness  of  the  boundary 
layer  where  concentration  polarization  takes  place. 

3.  Theoretical  Studies  of  Mass-Transfer  Rate.  In  1965,  Brian  developed  a theo- 
retical model  for  predicting  the  mass-transfer  rate  in  a parallel-plate,  reverse  osmosis 
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system.'’  Ilis  model  involved  a niimerieal  solution  of  a simplified  version  of  the  steady 
state  eiiuation  of  eontinuity  of  the  nontransferred  solute.  Brian’s  equation  is  as  fol- 
lows: 


3(U  • C)  ± 
dX  3Y 


D9C 

3Y 


0= 


(3) 


where:  U 

X 
C 

Y 

V 
1) 


axial  velocity 
axial  distance 

concentration  in  the  liquid  stream 
lateral  distance 
lateral  velocity 

effective  diffusivity  or  reciprocal  Peclet  Number  at  initial 
conditions 


The  terms  in  liquation  (3)  are  in  direct  correspondence  with  the  theory  of 

8(U  • C) 

concentration  polarization  discussed  earlier.  The  first  term, — — , represents  forced 

oX 

convection  of  solute  axially  through  the  cell.  The  first  part  of  the  second  term, 
3 

— : ( V • C),  indicates  bulk  convection  of  solute  toward  the  membrane;  and  the  second 
D3C 

part  ol  the  second  term.  . represents  molecular  diffusion  of  rejected  solute  away 
from  the  membrane  surface. 


In  related  studies,  Brian,®  Sherwood,  ft  a/.*  and  Gill,  el  al.’’  * developed 
slightly  different  models  modifying  the  boundary  conditions  and  solution  techniques 
of  Equation  (3).  Unfortunately,  none  of  the  studies  included  experimental  verification 
of  the  predicted  results. 


P.  L.  T.  Brian.  “Concentration  Polarization  in  a Reverse  Osmosis  Desabnation  with  Variable  I'lux  and  Incom 
plete  Salt  Rejection."  InJ.  F.nft.  Chem.  Fund.,  4.  439-45  ( 1 965). 

^ P.  L.  T.  Brian.  “Concentration  Polarization  in  a Reverse  Osmosis  System,"  Ind.  Fng.  Chem.  Fund.,  5,  148-149 
(1965). 

^ T.  K.  Sherwood,  P.  L.  T.  Brian.  R.  K.  Fisher,  and  L.  Dreshcr.  “SaH  Concentration  of  Phase  Boundaries  in  De* 
salination  by  Reverse  Osmosis,"  Ind.  Fug.  Chrm.  Fund.,  4,  1 1 3*18  ( 1965). 

^ W.  N.  (iill,  C.  Tien,  and  1).  W.  Zeh.  “Boundary  Layer  Lffects  in  Reverse  Osmosis  Desalination,”  Ind.  F.nft. 
Chrm.  Fund.,  5.  367-70  (1966). 

W.  N.  (lill,  C.  1 ien,  and  1).  W.  Zeh.  “Concentration  Polarization  Kffects  in  a Reverse  Osmosis  System  “ Ind. 
Fnf(.  Chvm.  Fund..  4.433-39  (1965). 


( 


5 


With  the  increasini’  interest  in  the  application  of  membrane-separation  proc- 
esses. several  theoretical  and  experimental  methods  of  redncint!  concentration  polari/a- 
lion  effects  have  been  studied.  Sherwood,  cl  nl.'*  performed  an  experimental  and  theo- 
i.'iic.*!  invc^^igation  of  a rotating,  cylindrical,  reverse  osmosis  apparatus  tor  desalination 
.md  i.ibtamed  partial  confirmation  of  their  previously  developed  theory,'®  Rotation  ot 
the  membrane  resulted  m total  llux  increases  of  ap)iroxmiately  lU  percent.  Zeh  and 
lull''  published  a theoretical  model  of  a rotating-disk  system,  implying  but  not 
tlLiantitatively  reporting  the  beneficial  effects  ol  rotation  on  the  mass-transt'er  rate 
Increa.sed  velocities  in  laminar  tlow  are  also  known  to  decrease  concentration  iiolari/a- 
tion  by  crr.iting  greater  shear  rates,  and  turbuient  How  adds  the  mechanism  ot  eddy 
mixing.  With  the  highly  permeable  membranes  available  today,  the  fluid  velocities  re- 
quired to  significantly  reduce  concentration  polarization  become  very  lar.ae.  and  energy 
consiimptiim  and  lluid  recirculation  rates  required  to  achieve  reasonable  recovery  be- 
come prohibitive.'* 


l ien  and  tiill'-’  theoretically  investigated  the  use  of  impermeable  “relaxation 
se  tioiis”  to  induce  axial  diffusion  of  the  non  transferred  solute,  but  their  predictions 
indicated  only  a small  increase  in  the  steady  state  tlux  obtainable  with  an  equivalent, 
entirely  penneable  area.  Huffman'''  studied  the  effect  of  natural  convection  on  con- 
centration polarization  by  using  a horizontal,  parallel-plate  ultralilter.  He  operated  a 
single  membrane,  first  in  the  lower  position  and  then  m the  upper  position,  and  ob- 
tained semitheoretical  correlations  for  his  forced-only  and  forced-plus-natural  convec- 
tion data.  The  study  demon.  ..^led  that  natural  convection  etfects  can  increase  the 
steady  state,  mass-transfer  rate  by  as  much  as  200  percent  over  comparable  conditions 
with  (orced  convection  only. 

In  1^72.  (jooding'''  experimentally  validated  his  dynamic  model  of  forced- 
convection  ultrafiltration  in  thin,  horizontal  channels  which  requires  a numerical 
solution  ot  tlie  unsteady  state  continuity  equation  for  the  nontransferred  solute.  He 

® T K.  .Slii-rvvood.  I’.  L,  T lirian,  and  R.  i;.  1 istu-r.  “Dcsaliiution  b\  Reverse  Osnuisis."  tnd.  h im.  C.hrm.  fund.. 
6.  2-i:  (1967). 

T K.  Sherwood.  P.  L.  t.  Brian.  R.  t . I isher.  and  L.  Diesher  “Self  Concentration  ol  Phave  Boundaries  in  He- 
vaiination  by  Reverse  Osmosis."  Ind.  h ('.hem.  i'und..  4.  113-18(1 965  i, 

* ^ I).  V^.  Zeh  and  W.  N.  CfiU.  “Convective  Diltusion  in  Rotating  Disk  Sysiems  with  an  Imperfect  Semi-permeable 
Interlace.”  1 mer.  In.st.  (.henu  I nfi.  jounuth  14.  715-19  (I968l. 

\ S.  Micliaeis.  "New  Separation  1 echnique  for  ihe  CPI.”  64.  No.  12.31-43(1968). 

C.  lien  and  W.  N.  Gill.  “The  Relaxation  of  ('oncentration  Polarization  in  a Reverse  Osmosis  Desalination  Sys- 
tem,” 1/ncr.  Inst  i'hett,.  t.rifi.  Jonrnal.  12.  722-27  (1966). 

W J.  Hutfinan.  “The  I ffcct  ot  l-orced  I nd  Natural  Convection  Dunne  the  l liraliliratiou  of  Protcin-Salme 
Solutions  in  I hln.  Horizontal  ('hanncls.”  I‘hh  DtssertotunL  CIcmson  I niversity  I ibrary.  ( lemson.  S (. . ( 1970). 

M C.  Gooding.  “A  Thcor»*lieal  and  1 xpcrimental  Comparison  of  Unsteady  and  Steady  State  lltratiltralion  in 
1 bin  Ilori/ontal  Channels.”  U.S,  Thesis,  CIcmson  University  ( 1972). 
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tlu.'n  iiscd  this  model  to  study  the  possibility  of  improving!  m.iss  traiisler  by  usiiip  some 
form  ot  unsteady  state  operation.  Ilipher  flu.xes  were  predieted  by  the  model  but  were 
not  experimentally  aehieved  in  his  investijration. 


4.  Practieal  Considerations  of  Permeate  Production.  I he  details  ol  the  deriva- 
tion of  the  praetieal  equation  of  concentration  polarization  control  of  ultraliltrate  lliix 
have  been  worked  out  by  others,'*  Simply,  the  relationship  is: 

f =■  K fn  (t  ^/Ch  I (4, 


where:  1'  = 

ultrafiltration  flux 

= 

concentration  of  the  gel  at  the  membrane  surlaee 

C’b  = 

concentration  ol  the  bulk  stream 

The  mass-transfer  coefficient  is  defined  m laminar  tlow  as: 

YO-33  |jU.67 
|,0..13  l".33 

(5i 

where:  V = 

velocity  along  the  membrane 

1)  = 

diffusivity  of  solvent 

1 

h 

height  of  the  channel  in  which  the  Hiiid  is  flowing 

L = 

length  of  the  channel. 

In  turbulent  How. 

the  mass-transfer  coefficient  is  defined  as: 

1 

75  j^fi.67 

I K = 

,,.i.;,s  ^(1.42 

where  v = the  kinematic  viscosity  of  the  I'luitl. 

While  these  eiiuations  are  useful  in  the  predication  of  effects  of  channel 
heipht  and  length  on  Htix,  irredietion  of  actual  system  operation  for  wastewater  treat- 
ment is  not  possible. 

S,  Mcssinjicr.  "l  ltraliltrali**n  .A  New  Proaws  ( i>ol.*’ 
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5.  Buckgruuml.  A matliematical  lunction  is  defined  as  a relationship  between  a 
set  of  X values  associated  witli  the  independent  variable,  x and  a set  of  y values  asso- 
ciated with  the  dependent  variable,  y.  Symbolically,  we  described  this  mathematical 
relationship  by; 

y = fix),  (7) 

The  function  can  also  be  displayed  by  a graph,  or  curve,  depicting  a set  of  data  [>airs 
(x,  y). 


When  we  are  given  a mathematical  equation  describing  a relationship  be- 
tween X and  y,  a curve  corresponding  to  this  relationship  can  be  easily  plotted.  On  the 
other  hand,  when  dealing  with  experimental  information,  the  opposite  is  true.  From  a 
given  set  of  points  (x,  y),  we  try  to  find  a function  which  generates  a curve  fitting  these 
data  points. 

For  the  purpose  of  this  study,  all  functions  are  assumed  to  be  continuous. 
This  means  that  for  every  value  of  x in  an  interval  (a,  b)  there  is  a unique  value  tor  y = 
fix)  and  that  a small  change  in  x produces  only  a small  change  in  y. 

6.  General  Curve  Fitting.  In  technological  applications,  measurements  are 
usually  made  of  one  quantity  while  a second  quantity  is  varied  with  the  results  being 
plotted  in  graphical  form.  The  problems  lie  in  not  being  able  to  make  measurements  at 
every  point  of  interest  and  the  possibility  of  measurement  error  in  the  measurements. 

To  overcome  these  problems  and  to  be  able  to  represent  the  results  in  a com- 
pact manner,  the  experimenter  often  tries  to  fit  the  experimental  measurements  to  a 
smooth  curve. 

The  form  of  the  curve  selected  will  of  course  depend  upon  the  type  of  mea- 
surements being  made  and  the  system  under  investigation.  A standard  function  often 
used  to  represent  a curve  is  a polynomial  function  ol  the  lorm. 

fix)  = ay  + a,  X -I- a2  x^  -I- . . . . a^,  x“.  (S) 

The  shape  of  the  resulting  curve  depends  both  upon  the  value  of  the  coefficients,  a. 
and  the  degree  of  the  polynomial,  u. 

7.  Least-Squares  Curve  Fit.  While  trying  to  fit  a curve  to  a set  ol  data  points,  it 
is  not  possible  to  have  the  curve  pass  through  each  point.  The  most  desirable  condition 
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is  the  curve  coming  as  close  as  possible  to  all  the  data  points.  If  the  set  of  data  points 

is  denoted  as  (x, , y,  , yj ) (x^,  y^ ).  where  is  the  value  of  the  independent 

variable  when  the  measurement  y^  was  made,  then  the  distance  of  the  i data  point  from 
the  curve  f(x)  is: 

= Vi  - t'(x,)  = y^  - (a„  + a,  Xj  + aj  Xj^  + . . . a^  x''.  (9) 

This  distance  represents  the  “error”  between  the  value  given  by  the  curve  and  the  mea- 
sured value. 

The  common  formula  that  measures  the  difference  between  the  curve  and 
the  data  is  the  sum  of  the  squared  errors.  That  is: 

C=  I y -f(x.)^  . (10) 

1-1  ‘ ‘ 

To  minimize  the  total  error,  the  coefficients  of  a through  a,^  are  selected  in  such  a 
manner  as  to  minimize  C.  This  function,  having  the  smallest  C value,  is  called  a least- 
squares  approximation.  The  conditions  that  must  be  satisfied  are  shown  as  follows 


for  u = 1,2,  and  3: 

Case  1 : u = 1 Linear  fit,  f (x)  = ao  -i-  a,  1 x:  (11) 

aov -t  a,  2xj  = Sy^  (12) 

ao  2 Xj  + a,  2Xj^  = 2 Xj  y^  (13) 

Case  2:  u = 2 Quadratic  fit,  fix)  = ao  + a.  y^  + a2  xT  ( 14) 

ao  + a,  2 Xj  + a22  xT  = 2 y.  (15) 

ao2  Xj  + a,  2 Xj^  + aj2  xT  = 2xj  y.  (Ki) 

ao2  Xj^  + a,  2 x^’  -(■  a22  xT  = 2xT  y.  (17) 

Case  3:  u = 3 Cubic  fit.  fix)  = ao  + a,  Xj  + a2  xT  -i-  a^  x^^  (IK) 

ao  + a,  2 Xj  + a22  xT  + aj2  Xj^  = 2y^  (19) 

ao2  Xj,  + a,  2 x3  + a,^  2 Xj^  + aj2  x^^  = 2x.  y^  (20) 
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I lu'SC 
C'lisi.'  1 ; 


+ a. 

v.x,^ 

+ 11:-  x," 

-r  aji  x,^  = lx,=  y, 

(21) 

+ a, 

vx; 

-)-a,,l  x,*’  = iix,^  y,. 

(22) 

in  turn,  can  be  solved  lor  the  difCerenl  coeft'icienls  of  a as  shown  in 


t ase 


y,)  (-X,  > 
) - (iix,-  I 


(23) 


a ml 


11 1 


V (ilx,  V|)  - (iix,)  (ly,) 

(YxTTTiTxTl 


( M) 


Or,  more  simply,  the  simultaneous  eciuations  can  be  solved  directly  by  the  method  ol 
elimination  or  by  the  use  ol  determinants. 


It  can  readily  be  seen  that  the  higher  the  value  of  u the  better  the  probable 
curve  fit  for  nonlinear  data.  It  can  also  be  seen  that  the  hipher  the  value  ol  u the  more 
difficult  it  will  be  to  solve  for  the  eoeftieients  ol  a.  Theretore.  other  limctions  are 
tried  and  compared  with  results  from  polynomial  fits. 

One  function  considered  is  an  exponential  1 unction  as  lollows: 

V = ae^"  . 


,Vt  first,  this  appears  to  be  a difficult  function  to  work  with;  however,  by  takine  the 
Vn  of  both  sides  its  form  becomes  a linear  function  similar  to  Case  1 as  shown' 

Vn  Y = tn  a + b X . (261 

H.  Correlation.  Once  a curve  or  a line  is  fitted  to  the  data,  the  next  problem  is 
to  determine  how  well  the  data  is  actually  fitted  by  the  curve.  .‘\  lair  idea  is  gained  b\ 
inspection  of  the  curve  and  the  data  points;  however,  usually  a more  objective  evalua- 
tion is  re(|uired.  In  studying  the  data,  an  as.se.ssmejit  must  be  itiade  ol  the  probability 
of  random  variables  rather  than  the  independent  variable  contributing  to  the  variabilitv 
of  the  data. 

b.ssentially . the  problem  is  analysis  ol  variance;  ami  the  (luantity  to  be  anal- 
y/ed.  the  total  sum  of  the  squares,  is  given  by: 

1 (y|  - yi^  (27) 

r 1 
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(28) 


v\liere:  y = the  mean  of  the  y values  jiiven  by  y = iy 

1=1  ‘ 


II'  the  iruiependent  variable  were  the  only  thing  that  alfeeted  the  dependent 
variable,  the  data  points  would  all  fall  on  a straight  line  assuming  a linear  relationship. 
Since  this  is  rarely  true,  other  factors  do  tend  to  affect  the  dependent  variable.  'I'he 
variation  caused  by  these  other  factors  is  usually  measured  by  the  sum  of  the  stjuares 
ol  the  vertical  deviations  from  the  point  to  the  line  given  by: 


-vl’  C9| 

1 = 1 


where  y = the  individual,  predicted  y values. 


riiis  is  itleiuical  to  eiiuation  (10)  or  f(Xj).  The  difference  betweeen  eipia- 
tioii.  (27)  and  (29)  measures  the  variation  of  the  y’s  that  can  be  attributed  to  x.  This 
valiK  divided  by  equation  ( 27 ) gives  the  proportion  of  the  total  value  of  the  y’s.  or  de- 
pendent variable,  that  can  be  attributed  to  the  relationship  with  x and  is  called  the  co- 
efficient of  determination.  The  square  root  is  called  the  coefficient  of  correlation  with 
the  sign  of  this  being  the  same  as  that  of  the  coefficient,  b.  m the  equation  of  the  least- 
squares  line.  These  eciuations  are  as  follows; 


/ 

.2  _ -(Vi* - y>^  "“<^1  ■ • 


(30) 


and 


, lAoi-y)-^ 
r 2;(y, -y)2 


* -2;(y,  -y)' 


(31 ) 


where-*  = the  coefficient  of  determination 
r = the  coefficient  of  correlation. 

The  value  of  r is  the  most  wiilely  used  measure  of  the  strength  of  the  linear 
relationship  between  two  variables.  It  indicates  the  goodness  of  fit  of  a line  or  curve 
fitted  by  the  method  of  least  squares  and  this,  in  turn,  tells  whether  or  not  it  is  reason- 
able to  say  that  a linear  relationship  (correlation)  exists  between  the  dependent  and  in- 
dependent variables. 
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IV.  I^XPliRIMl  NTAL 


Data  for  tliis  study  came  from  a Department  of  rransportation  (DOT)  report 
dealing  with  the  renovation  of  shower  wastewater.'^  The  purpose  of  the  DOT  report 
was  to  evaluate  the  effectiveness  of  ultrafiltration  for  the  renovation  of  shower  waste- 
water.  Different  membranes  and  membrane  configurations  were  utilized  on  the  various 
treatment  systems. 

9.  Membrane  Configuration.  The  current  study  evaluated  two  different  sys- 
tems using  two  different  membranes  in  two  different  membrane  configurations.  These 
membrane  configurations  were  hollow-fiber  and  spiral-wound.  Both  configurations 
employ  minimal  channel  height  which  according  to  equations  (5)  and  (6)  increases  the 
mas.s-transfer  coefficient. 

The  hollow-fiber  membrane  module  contains  660  fibers  of  0.045-inch  bore 
diameter.  Hach  fiber  consists  of  a copolymer  support  with  the  noncellulosic  membrane 
approximately  10.0  microns  in  thickness  coated  on  the  inside  of  this  support  (l  igure 
4).  The  fibers  are  encapsulated  at  both  ends  in  a 3-inch-diameter  by  25-inch-long, 
plastic  cartridge.  An  end  view  of  the  cartridge  is  shown  in  Figure  5.  The  total  mem- 
brane area  is  15  scjuare  feet.  The  feed  stream  of  raw  wastewater  flows  through  the 
inside  of  the  fibers,  and  the  permeate  collects  in  the  shell. 

The  spiral-wound  membrane  consists  of  a sheet  22  inches  wide  with  a 
support  backing  and  0.030-inch-thick  vexar  spacer  material.  This  entire  system  is 
then  rolled  in  a jelly  roll  fashion  around  a permeate  collection  tube.  Wastewater 
passes  through  the  channel  made  by  the  vexar  spacer  material  which  also  serves  as  a 
turbulence  promoter.  The  penneate  passes  through  the  membrane  and  is  channeled 
by  the  support  material  to  the  penneate  collection  tube.  Wastewater,  called  waste 
brine,  exits  the  end  of  the  channel  made  by  the  vexar  spacer  material.  The  module 
contains  approximately  30  square  feet  of  membrane  area.  Fnd  and  side  views  of  the 
module  are  shown  in  Figure  b. 

10.  System  Description.  A simplified  fiow  diagram  of  the  hollow-fiber  system 
operated  in  batch  mode  is  shown  in  Figure  7.  The  wastewater  is  pumped  from  a feed 
tank  through  a bag  prefilter,  fhe  ultrafiltration  system  is  operated  at  a pressure  of 
approximately  15  psig  and  an  inlluent  fiow  rate  of  16  gallons  per  minute  at  a velocity 
of  6 feet  per  second.  While  permeate  is  collected  in  the  shell  of  the  module,  the  con- 
centrate is  recirculated  to  the  feed  tank.  The  system  also  had  the  capability  of  being 


I)  S.  Lent.  “Renovation  ot  Waste  Shower  Water  by  Membrane  I ilfration.”  Interim  Report.  I'.S.  Department  ol 
Transportation  (19771. 
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Figure  7.  Flow  diagram  of  the  hollow-fiher  system. 
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uliidi  allows  tla'  iK'niicak-  to  lu'  puiii|ii.'i.l  I'roin  the  outside  oT  IIk’  I'ihei  to 
the  inside,  thus  displaeing  any  itel  layer  buildup  on  the  ineinbrane  surtaee. 

Daily  operation  eonsisted  of  (i  to  8 hours  of  eonliifious  i p.eration.  lollow- 
int;  the  eiul  of  eaeh  day’s  run.  the  system  was  baekl'lushed  with  about  20  gallons  of 
permeate.  On  the  followini:  rlay,  the  direetion  of  flow  throuiih  the  eartridee  was  re- 
versed, thus  eliminatinu  any  buildup  of  partieulate  material  at  tl"’  head  or  entranee,  of 
the  membrane  module.  The  system  eoneentrates  a tank  ofwa'  to  about  15  percent 
of  its  original  volume.  I'he  tank  is  then  dum|ied  aiul  refilled  with  fresh  wastewater. 

A simplified  How  diagram  of  the  system  with  the  spiral-wound  membrane 
configuration  operated  in  feed  and  bleed  mode  is  shown  in  f igure  8.  Wastewater  is 
intmped  through  two  prefillers  into  the  eiretdation  [nimp  which  maintains  the  velocity 
in  the  circulation  loop.  Within  the  circulation  loop,  two  membrane  cartridges  are  con- 
nected in  series.  The  flow  rate  in  the  loop  is  10  to  15  gallons  per  minute.  Concentrate 
from  within  the  circulation  loo|i  is  purged  via  a solenoid  valve  operated  periodically  to 
maintain  the  desired  tiegree  of  concentration  in  the  loop. 

This  system  was  also  operated  b to  8 hours  daily.  I'he  degree  of  concentra- 
tion maintained  in  the  loop  is  appro.xiinately  98  percent. 

I 1.  Data  Analysis.  Two  hundred  and  twenty-two  hours  of  operati  /iial  data  for 
the  hollow-nhcr  system  and  1 18  hours  of  data  for  liie  .sp>ira)-vvound  system  were  anal- 
yzed. Data  analysis  was  then  performed  on  the  first  118  hours  of  data  from  the 
hollow-fiber  system  to  give  an  ecjual  time  base  for  system  comparison.  The  data  for 
the  hollow-fiber  system  was  then  analyzed  on  a daily  basis  because  at  the  end  of  each 
day's  run  the  system  was  backnushed  or  cleaned  Cu'iving  a significant  increase  in  per- 
meate production. 

The  data  analysis  consists  of  first,  second,  and  third-order  polynomial  curve 
fits  and  an  exponential  curve  fit.  The  analysis  was  performed  on  a Control  Data  Cor- 
I poration  6600  Computer  using  a Tektronix  4014  Clraphics  Terminal,  The  analy  sis  con- 

sists of  a scatter  gram  plot  of  the  data,  a curve  generated  to  best  fit  the  data,  the  epua- 
tion  of  the  curve,  and  the  st|uare  of  the  multiple-correlation  coefficient. 

V.  DLSCUSSION 

12.  Analysis  of  Results.  I xperimenlal  values  used  in  this  sliiily  are  shown  in 
Appendix  A.  ’I’he  computer  program  useil  for  the  itata  analysis  is  shown  in  .Appendix 
Ik  Table  I gives  the  coefficients  of  determination  for  various  orders  of  fit.  fable  2 
gives  the  coefficients  of  determination  for  the  time  periods  between  backnushes  for 
the  hollow-fiber  system. 
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Figure  8,  Flow  diagram  of  the  spiral-wound  system 


1 ;ible  I . ( ootTicicMits  of  Determination  lor  Various  Orders  ol'  l it 


System 

Orde 

r of  Fit 

1st 

2nd 

3rd 

L 

1 lollow-Fiber, 

.02(1 

.201 

.302 

.015 

All  (222  Hr) 

Hollow-Fiber, 

.3.50 

.514 

.732 

.415 

Part  ( 1 18  Hr) 

Vexar  Spiral- 

.501 

.837 

.944 

.(■>11 

Wound  (118  Hr) 

Results  of  the  runs  performed  on  tlie  total  data  for  both  systems  are  shown 
in  Appendix  C along  witli  the  results  for  1 18  hours  of  operation  on  the  hollow-fiber 
system. 

Appendix  D contains  the  daily,  segmented,  or  between-backnushing,  data 
for  the  hollow-fiber  system.  Appendix  K contains  segmented,  fitted  data  for  the  spiral- 
wound  system. 

Table  I shows  that  the  best  curve  fits  are  obtained  with  third-order  poly- 
nomials. Thus,  the  general,  best-fit  equation  is  of  the  form  Y = a x^  + b + c x + d. 
This  is  further  demonstrated  in  Table  2 where  once  again  third-order  polynomials  give 
the  best  fit.  The  trend  for  best  fits  shows  that  second-order  fits  follow  in  fitting  the 
data,  followed  by  exponential  fits,  followed  by  a first-order,  or  linear,  fit, 

fxamination  of  the  graphs  for  the  total  data  in  Appendix  C.  however,  shows 
that  the  curves  generated  by  second-  and  third-order  polynomials  describe  unrealistic 
results.  The  second-order  fits  in  all  cases  show  a maximum,  a minimum,  and  then  a 
maximum.  This  contradicts  the  phenomenon  of  concentration  polari/.ation  described 
earlier  in  that  the  boundary-layer  accumulation  would  have  to  increase  to  a point 
(causing  a decrease  in  tlux)  then  decrease  at  approximately  the  same  rate  (causing  an 
increase  in  flux). 

The  third-order  fits  go  one  step  further  from  the  maximum/minimum/maxi- 
mum and  add  a minimum.  It  is  possible  for  the  iTux  to  decrease  and  increase;  how- 
ever, the  rate  of  increase  is  always  rapid  and  is  achieved  when  the  boundary  gel  layer  is 
suddenly  removed  either  naturally  or  by  cleaning. 

The  empirical  model,  therefore,  would  be  an  equation  of  hign  order:  thus, 
by  having  the  iiredicted  equation  oscillate  above  and  below  the  data  enough  times,  the 
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Tahk'  2.  ('oc('fideiit.s  of  Detcrminulion  lor  1 iiTie  I'eriods 
Ik-twfcn  Backnusliings  tor  the  Hollow-l-'ibcr  System 


S -gment 

Order  of  Fit 

1st 

2nd 

3rd 

L 

0 to  8.b 

.720 

.927 

.972 

.827 

8.6  to  15.4 

.830 

.983 

.989 

.861 

15,4  to  19.6 

.796 

.923 

.925 

.792 

20.2  to  44.2 

.870 

.917 

.958 

.901 

44.9  to  48.9 

.453 

.555 

.599 

.456 

52.2  to  60.0 

.779 

.926 

.932 

.782 

59.9  to  90.7 

.867 

.923 

.973 

.910 

91.4  to  95.3 

.657 

.849 

.915 

.710 

98.4  to  100.3 

.997 

1.000 

1.000 

1.000 

101.5  to  104.9 

.808 

.946 

.994 

.863 

106.1  to  1 14.1 

.480 

.573 

.584 

.555 

115.6  to  122.9 

.004 

.020 

.234 

.007 

124.0  to  130.1 

.697 

.876 

.880 

.684 

1 30.6  to  1 36.9 

.894 

.978 

.978 

,908 

137.9  to  145.1 

.781 

.943 

.945 

.785 

145.1  to  150.4 

.885 

.887 

.931 

.898 

150.6  to  157.3 

.654 

.896 

,929 

.682 

157.7  to  164.8 

.291 

.782 

.896 

.286 

164.8  to  170.5 

.456 

.506 

.593 

.455 

171.3  to  176.1 

.802 

.843 

.971 

.815 

176.3  to  184.5 

.954 

.967 

.971 

.969 

185.3  to  190.9 

.732 

.734 

.895 

.721 

192.3  to  199.6 

.580 

,940 

.995 

.585 

200.2  to  207.6 

.213 

.213 

.834 

.207 

208.1  to  215.7 

.539 

.769 

.831 

.543 

216.5  to  222.4 

.840 

.851 

.923 

.723 

equation  very  nearly  fits  tlie  line.  However,  tlte  problems  with  this  would  be  the 
large  number  of  eonstants  needed  to  define  the  equation  and  the  computation  of  the 
constants.  Therefore,  without  the  data,  prediction  of  any  results  would  be  purel.v 
arbitrary. 

I.T  Proposition  of  Theory.  Consideration  of  the  actual  system  dictate.s  that  ini- 
tially or  early  in  the  run  the  system  is  at  unsteady  state.  This  is  true  because  at  time 
zero  there  is  no  boundary  layer  present  which  is  evident  by  observing  the  initial,  rapid 
rate  of  fiu.s  decline  which  is  pro[iortional  to  the  rate  of  buildup  of  boundary  layer. 


i 


:o 


I'liis  buildup  is  dircciK'  pro|n)rtioiKil  to  tlu.'  solute  concentration  because  almost  all  ol 
the  solute  removed  from  the  waste  stream  is  I'ormme  the  boundary  layer. 


t 
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As  the  boundary  layer  develops,  the  transport  of  material  Irom  the  boumlary 
layer  to  the  bulk  solution  begins  to  be  the  controlling  lactor  lor  iiermeate  production, 
rite  amount  of  this  lateral  diffusion  is  more  proiiortional  to  the  system's  design  and 
operational  parameters  than  to  the  solute  concentration.  1 his.  ol  course,  is  true  only 
under  the  assumption  that  permeate  production  is  proportional  to  the  amount  ol 
boundary  layer. 

14.  Proposition  of  the  Model.  Under  this  [iremise.  a two-part  model  was  tried. 

1 he  first  part  would  cover  the  unsteady  state  condition  ot  rapid  llu.\  decline  and  the 
second  part  would  cover  the  condition  where  dillusion  away  (rom  the  membrane  is 
the  controlling  factor.  A more  detailed  approach  would  also  consider  the  transition 
period  where  permeate  production  is  a function  ol  both  parts,  ilowever,  since  specilic 
boundaries  for  this  transition  state  would  be  very  hard  to  establish,  the  two-|iart  model 
was  chosen. 

The  general  procedure  for  using  the  model  is  to  visualK  inspect  the  complete 
plots  of  the  data  and  the  data  itsell.  Then,  the  data  is  separated  into  two  portions 
which  tentatively  satisfy  the  conditions  of  the  model  and  perform  curve  fits  on  both 
portions  of  the  data. 

15.  Application  of  the  .Model  to  the  Hollow-Fiber  Syslcm.  laxaminalion  of  the 
plots  and  data  of  the  hollow-fiber  system  in  Appendices  A and  C tails  to  reveal  the  two 
parts  described  earlier.  What  the  examination  does  reveal,  however,  are  many  seg- 
ments. each  exhibiting  a rapid  rate  of  flux  decline  with  varying  degrees  ot  Hux  recovery 
between  each  segment.  Therefore,  it  appears  that  this  system's  data  does  not  fit  the 
model. 

However,  when  the  daily  opeiation  of  the  system  is  considered,  it  is  clear 
that  the  data  does  fit  the  model.  This  is  true  because  each  day  the  system  was  back- 
tlushed.  This  operation  consisted  of  pumping  permeate  backward  through  the  mem- 
brane thus  removing  the  boundary  layer.  Since  the  boundary  layer  was  removed  each 
day,  it  never  had  a chance  to  build  up  enough  to  become  the  controlling  factor  in  per- 
meate production.  Therefore,  the  system  was  always  in  an  unsteady  state  condition, 
and  only  the  first  portion  of  the  model  should  be  applied  to  the  system's  data.  Verifi- 
cation of  the  model  consisted  of  segmenting  the  data  into  time  periods  between  back- 
tlushes  and  then  performing  curve  lits.  The  coetficients  ol  determination  tor  the  litted 
curves  are  shown  in  Table  2,  and  the  actual,  litted  curves  are  graphically  shown  in 
Appendix  D. 


A comparison  of  tlic  codTicicnts  of  determination  from  the  entire  run  and 
tlie  first  I 18  hours  of  the  run  witli  the  coefficients  from  the  seftmented  data  sliows 
mucli  better  correlation  witli  the  segmented  data.  Review  of  tlie  second-  and  tliird- 
order  plots  in  Appendix  1)  also  shows  these  plots  to  be  somewhat  unrealistic  because 
of  their  cyclic  nature.  I'irst-order  fits  can  also  be  eliminated  since  the  trend  in  all  cases 
is  r -ver  linear.  I'herefore.  the  exponential  fit  best  fits  the  segmented  data  for  the 
hollow-fiber  system.  The  coefficients  as  well  as  the  coefficients  of  determination  for 
the  various  segmented  sections  are  giveti  in  Table  3. 


fable  3.  Txponential-Curve-Tit  Coefficients  and  Coefficients  of 
Determination  for  the  llollow-F'iber  System 


Segment 

a 

b 

r' 

0 



8.6 

79.0 

-0. 1 04 

.827 

8.6 

- 

15.4 

102 

-0.070 

.861 

15.4 

- 

19.6 

172 

-0.098 

.79  2 

20.2 

- 

44.2 

43.9 

-0.024 

.901 

44.9 

- 

48.9 

60.6 

-0.022 

.456 

52.2 

— 

60.0 

56.0 

-0.018 

.782 

59.9 

90.7 

51.1 

-0.014 

.910 

91.4 

— 

95.3 

1.35  x lO'' 

-0.063 

.710 

98.4 

— 

100.3 

2.47  X 10'  ' 

-0.230 

1 .000 

101.5 

- 

104.9 

1.42  X 10“ 

-0.151 

.863 

1 06. 1 

— 

1 14.1 

3.15  X lO'* 

-0.069 

.555 

1 15.6 

- 

122.9 

18.1 

-0.003 

.007 

124.0 

- 

130.1 

1.61  X 10'’ 

-0. 1 09 

.684 

130.6 

— 

136.9 

1.74  X 10“ 

-0.050 

.908 

137.9 

— 

145.1 

3.28  X 10^ 

-0.035 

.785 

145.1 

- 

1 50.4 

4.40  X 10“ 

-0.201 

.898 

150.6 

_ 

157.3 

3.87  X lO’ 

-0.0932 

.682 

157.7 

— 

164.8 

174 

-0.01  15 

.286 

1 64.8 

— 

1 70.5 

169 

-0.0212 

.455 

171.3 

176.1 

1.28  X 10^ 

-0.0225 

.815 

1 76.3 

— 

184.5 

3.09  X 10“ 

-0.0504 

.969 

185.3 

190.9 

7.12  X 10" 

-0.0800 

.721 

192.3 

- 

199.6 

1.67  X 10-" 

-0.0936 

.585 

200.2 

- 

207.6 

417 

-0.0124 

.207 

208.1 

215.7 

4. 1 1 X 1 0“ 

-0.0341 

.543 

216.5 

— 

222  4 

1.85  X 10'  ’ 

-0.165 

.723 

where: 

y = 

ae*"  and  r^ 

= coefficient  of  determination. 
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16.  Application  of  the  Model  to  the  Spiral-Wound  System.  I xainination  of  the 
graphs  in  Appciuli.x  ('  for  this  system  reveals  an  initial,  rapid  rate  of  llux  decline  h)l- 
lowed  by  an  almost  linear  leveling  out  of  the  data.  Therefore,  this  data  a[.pears  to  lit 
the  proposed  theory  rather  well.  To  test  this,  the  data  was  segmented  into  two  sec- 
tions, and  curve  fits  were  performed  on  each  section.  The  first  segment  was  531,9  to 
555,3  and  the  second  segment,  the  remainder  of  the  data,  556,2  to  649.6,  This  was 
called  (iroup  1 on  which  first-order  and  exponential  curve  fits  were  ruti.  This  process 
was  continued  for  Groups  2 and  3,  each  of  which  segmented  the  data  in  a dilTerent 
place.  These  curve  fits  are  shown  in  Appendix  12  The  resulting  coefficients  and  coef- 
ficients of  determination  for  the  fits  are  shown  in  Table  4. 

Table  4,  Tirst-Order  and  Txponential-Gurvc-T'it  Coefficients  and 
Coefficients  of  Determination  for  the  Spiral-Wound  Sysfem 


Group 

No. 

Segment 

1 St  Order 

Fxponemial 

a 

b 

r2 

*a 

b 

r2 

1 

531 .9 

555.3 

-.13 

76.8 

.910 

3.4  X lO'" 

-.060 

.952 

T 

531.9 

560.2 

-.12 

66.2 

.900 

6.8  X 10' 2 

-0.057 

.963 

3 

531.9 

- 575.0 

-7.3  X 10-2 

42.3 

.800 

2.5  X 10'° 

-0.43 

.909 

1 

556.2 

- 649.6 

-4.52  X 10'2 

3.29 

.543 

38.0 

-0.007 

.535 

2 

561.1 

- 649.6 

-3.70  X I0'2 

2.77 

.480 

19,7 

-0.006 

.463 

3 

576.0 

- 649.6 

-6.30  X 10-“ 

.864 

.067 

1.02 

-0.001 

.057 

Where : y = ax  + b for  first  order 

y = ae**^  for  exponential 

r^  = coefficient  of  determination  t 

* calculated  using  Hewlett  Packard  model  65  calculator  with  StatPack  lixponential  I 
Curve  Fit  Program 


Group  2 data  for  the  initial,  or  unsteady  state,  gives  the  best  fit.  It  is  no 
surprise  that  the  best  fit  comes  from  the  exponential  equation  with  a rather  high  co- 
efficient of  determination.  As  for  the  steady  state,  or  latter  portion,  of  the  data,  the 
optimum  tends  to  be  closer  to  Group  1.  In  this  ease,  the  first-order  fit  gives  a better 
fit  than  the  exponential  fit.  Examination  of  the  daily  operation  data  shows  essential- 
ly constant  daily  permeate  production;  however,  on  different  days  this  figure  fluc- 
tuates. The  coefficient  of  determination  is  low  and  appears  to?;ct  worse  as  less  data  is 
considered. 

The  reason  for  this  can  be  seen  by  comparing  the  curve  fits  for  the  total  data 
found  in  Appendix  C and  the  curve  fits  of  the  segmented  data.  Visual  examination  of 
the  560-  to  650-hour  range  in  the  graph  of  the  total  data  on  this  system  reveals  an  al- 
most linear  area  of  the  curve.  Compare  this  with  the  561.1-  to  649.6-hour  segmented 
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graph  round  in  Appendix  H.  Since  tlie  permeate  production  is  mueli  less  variable  here 
than  for  the  whole  system,  the  y axis  is  greatly  expanded  here  as  compared  to  the 
whole  system,  riuis,  the  amount  of  error  is  expanded  as  less  data  is  considered  not 
only  in  the  visual  presentation  but  also  as  witnessed  in  the  low  value  for  the  coefficient 
of  determination. 

The  apparent  discrepancy  between  Group  2 data  best  fitting  the  unsteady 
state  condition  and  Group  1 data  best  titting  the  steady  state  condition  can  be  partly 
attributed  to  Group  I’s  steady  state  data  being  more  voluminous  and,  therefore,  having 
a better  coefficient  of  determination.  However,  the  main  reason  for  the  discrepancy 
could  be  the  definition  of  the  boundary  conditions  for  steady  and  unsteady  states.  In 
actuality,  the  boundary  conditions  probably  lie  between  the  boundary  conditions  de- 
fined as  Group  1 and  those  defined  as  Group  2;  and,  second,  it  is  doubtful  that  there 
is  a discrete  boundary  between  the  two  states  resulting  in  a third,  or  transition,  state 
as  pointed  out  earlier. 

VI.  CONCLUSIONS 

17.  Conclusions.  Based  on  the  data,  this  study  concludes  that: 

a.  A single,  empirical  equation  cannot  describe  permeate  production  by 
idtrafiltration  of  shower  wastewater. 

b.  A model  based  on  unsteady  and  steady  state  phenomena  of  concentra- 
tion polarization  can  empirically  describe  the  ultrafiltration  process. 

c.  The  hollow-fiber  system  can  be  described  by  a single,  empirical  equa- 
tion of  exponential  order  because  the  system  as  studied  was  never  allowed  to  reach 
steady  state.  Therefore,  the  unsteady  state  equation  described  the  system’s  operation. 


d.  The  spiral-wound  membrane  system  can  be  described  empirically  by  an 
exponential  equation  during  the  boundary  layer  formation  and  by  a linear  equation 


during  steady  state  operation. 
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DATA 


riic  following  pages  are  a eoniputer  listing  of  the  data.  Line  numbers  100  through 
2010  are  the  data  from  the  hollow-fiber  system.  The  first  number  is  the  line  number 
and  is  followed  by  an  equal  sign.  The  second  number  is  the  time  factor  expressed  in 
hours,  and  the  third  number  is  the  llux  expressed  in  gallons  per  square  feet  of  mem- 
brane area  per  day.  After  Line  2010,  the  line  numbers  begin  again  at  100  and  continue 
through  1380,  This  is  the  data  from  the  spiral-wound  system.  Time,  again,  is  ex- 
pressed in  hours;  however,  it  is  not  normalized  and,  therefore,  begins  at  531.9.  Flux, 
this  time,  is  expressed  in  gallons  per  minute. 
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COMI’UTLK  l'KO(.K  \M 


An  intcraclivo  computer  (iropram  wa;,  wnllcn  in  I'orlran  IV  language  for  pcrl'orming 
tlic  curve  tits.  The  subroutines  for  the  graphics  portion  of  the  program  are  from  tile 
Tektronix  Advanced  (Iraphics  Package.  I'hese  include  IMTT,  Ti  RM.  RINTIT.  MMX, 
BI  LL,  SCURSR,  DLIM.X,  DLLMV,  1 RASl..  LINT.  SYMBl..  ( 111  ( R.  DSPLAY. 
( I’LOT.  MOVABS,  IIS  TRIN.  VLABl  L,  AMMODT.  and  I INI  11.  The  program  was 
run  on  a Control  Data  Corporation  (lOOO  Computer  using  a Tektronix  4014  (Iraphics 
Terminal. 

The  data  was  stored  on  a permanent  disk  file.  The  context  editor  was  used  to  segment 
the  data.  This  data  was  then  saved  as  TAPT,  1 without  line  numbers.  TAPli  1 was  re- 
wound and  the  compiled  version  of  the  program  was  executed.  Alternatively,  data 
could  be  entered  directly  from  the  terminal  by  connecting  TAPT:  1 atid  entering  the 
data  when  directed  by  the  eom|iuter. 

Upon  execution,  a scatter  gram  of  /.ero-order  fit  is  plotted.  Options  of  T,,  R.  L,  or  .N 
can  then  be  entered  as  explained  in  line  numbers  580  through  (ilO  for  the  graphical 
representation  desired.  The  output  consists  ot  the  order  of  fit,  graphical  data  presenta- 
tion, the  titled  curve,  the  cquution  of  the  curve,  and  the  coefficient  of  determination 
expressed  as  SQ  MULT  CORR  COTT-'T  IS.  In  some  cases,  a coefficient  of  the  equation 
is  given  hy  ********  indicating  insufficient  storage  capacity  for  the  size  number 
stored.  The  order  of  fit  is  not  shown  on  the  graphs  because  of  margin  limitations. 
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APPENDIX  C 


CURVE  FITS  FOR  SYSTEM  OPERATION 


The  lirst  tour  graplis  represent  the  total  data  for  tlie  liollow-fiber  membrane  system. 
The  next  four  graphs  represent  tlie  first  1 18  hours  of  operation  of  the  hollow-fiber  sys- 
tem. These  are  ineluded  for  comparison  of  the  two  systems  on  an  equal-time  basis. 
Tlie  last  tour  graphs  represent  the  total  data  for  the  spiral-wound  membrane  system. 
The  specific  order  of  the  graphs  is  first-,  second-,  and  third-order  followed  by  the  ex- 
ponential curve  fit. 
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CURVE  FITS  FOR  DAILY  OPERATION  OF  THE 
HOLLOW-FIBER  SYSTEM 


Tlie  graphs  in  this  section  represent  tlie  segments  sliown  in  Table  2.  liacli  segment  rep- 
resents tile  time  period  between  backnushes.  The  usual  length  of  time  covered  is  5 to 
8 hours.  However,  in  some  cases,  the  system  was  allowed  to  continue  running  through 
the  night.  The  length  of  these  runs  could  be  from  24  to  32  hours  depending  on  what 
time  the  following  day  the  system  was  backtlushed.  A total  of  26  segments  or  graphs 
is  considered  for  each  fit.  The  first  26  graphs  represent  first-order  fits  for  the  seg- 
ments. These  are  followed  in  turn  by  second-order,  third-order,  and  finally,  the  expo- 
nential fits. 
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APPINDIX  1 

tUUVl  MTS  TOR  SIOMTNTLI).  SPIRAL-WOUND  DATA 

11k.'  data  lor  the  spiral-wound  system  was  segmented  at  what  appeared  to  he  the  point 
where  the  system  reaehed  steady  state.  Tliree  different  points  were  pieked  whieh  re- 
sulted in  three  groups  of  data  eaeh  eontaining  two  segments.  The  points  ol  division 
are  555.3,  5()0.2.  and  575.  liaeh  group  contains  four  graphs,  a first-order  and  e.xpo- 
nential  fit  for  tiie  presumed  unsteady  state  portion  of  tlie  data  and  the  same  two  tits 
for  the  steady  state  portion  of  the  data.  It  sliould  he  noted  that  in  the  two  plots  ol 
the  steady  state  for  (Iroup  2 there  is  an  erroneous  data  point  at' 5h0  hours.  A quirk 
in  the  graphics  packages  made  iitclusion  of  this  point  necessary,  and  it  should  have  only 
a negligible  effect  on  the  etiuation  and  the  coefficient  of  determination. 
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